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1,  INTRODUCTION 

~^When  an  adaptive  array  is  implemented  digitally,  the  sample  covariance 
matrix  algorithm  provides  a  direct  method  of  computing  the  adaptive  weights 
and  rapid  convergence  independent  of  the  eigenvalues  of  the  covariance 

Dpt" 

matrix.  Previous  analyses  of  this  algorithm1-  J  have  assumed  that  the 
weights  are  computed  using  one  set  of  array  element  outputs  and  these 
weights  are  applied  to  later  array  outputs.  This  report  considers  the 
case  where  the  adaptive  weights  are  tested  against  the  same  set  of  data 
used  in  the  weight  computation. 

For  many  applications,  the  multiple  channel  sidelobe  canceller  is 
the  preferred  adaptive  configuration.  It  can  be  shown  that  the  sidelobe 
canceller  is  a  special  case  of  the  more  general  adaptive  array.^  In 
the  next  section  it  is  shown  that  the  general  adaptive  array  problem  can  be 
transformed  to  an  equivalent  sidelobe  canceller  problem,  a  form  which  is 
more  convenient  for  some  analyses.  It  is  also  shown  that  the  array  per¬ 
formance  is  independent  of  this  transformation,  and  that  the  effective 
weights,  output  S/N,  etc.,  can  be  computed  in  any  convenient  coordinate 
system  provided  the  transformation  of  coordinates  is  non-singular.  , 


2.  COORDINATE  TRANSFORMATIONS 

Let  X  denote  the  column  vector  of  array  element  outputs, 

XT  =  >x2’' ’ • ,xrP’  and  Sx  denote  the  corresponding  signal  vector.  The 

noise  covariance  matrix  of  the  array  outputs  is  try 
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Mx  =  E  X  X 


0) 


where  Mx  is  a  NxN  Hermitian  matrix  for  an  N  element  array,  E  denotes  the 
expectation,  +  the  complex  transpose,  and  all  noise  components  (but  no 
signal)  are  included  in  Mx- 


The  weights  which  maximize  the  S/N  ratio  are 


«x  *  "S’  sx 


and  the  corresponding  array  output  is 


t  t  -l 

l  =  Wy  x  =  s  M  ‘  X 

A  A  A 


With  optimum  weights,  the  output  S/N  ratio  is 


(2) 


(3) 


<wxV  t 

r  =  -*—■ — -  =  S..  M.. 1  S„ 
Wx  Mx  Wx 


XXX 


(4) 


Let  T  denote  any  non-singular  transformation,  and 


Y  =  T  X 


(5) 


In  the  new  coordinate  system, 


My  =  E  Y  Y+  =  T  Mx  Tf 


S  =  T  S 
y  x 


(6) 

(7) 


"y  •  m;'  sy 


(8) 
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Combining  5  through  8, 

Wy  =  (Tf)~ 1  Wx  (9) 

and 

Z  =  wj  Y  =  Wx  X  (10) 

i.e.,  the  output  of  the  array  is  unchanged  by  the  transformation. 

The  sample  covariance  matrix  in  X  coordinates  is 

K 

"x  ■  f  \  xk  Cl) 

X,  l_ 

where  Xk  denotes  the  ktn  independent  sample  of  array  element  outputs  and  K 
is  the  number  of  samples  in  the  estimator  of  Mx<  The  weights  based  on  a 
sample  covariance  matrix  are 

«x  ■ sx  cz) 

Replacing  Mx  and  M  with  the  corresponding  esimators,  M  and  M  ,  and  fol- 
lowing  the  analysis  of  the  preceding  paragraph,  it  can  easily  be  shown 
the  array  output  with  weights  based  on  a  sample  covariance  matrix  is  also 
independent  of  the  transformation  T.  Hence,  the  analysis  of  adaptive  array 
performance,  including  the  sample  covariance  matrix  algorithm,  can  be  per¬ 
formed  in  any  convenient  coordinate  system  provided  the  required  trans¬ 
formation  is  non-singular. 


/ 
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3,  PROBABILITY  DISTRIBUTION  OF  S/N 

For  any  arbitrary  adaptive  array,  the  input  vector  X  can  be  trans¬ 
formed  to  a  new  set  of  coordinates  in  which  the  signal  is  present  in  only 

one  component  and  the  noise  covariance  matrix  is  diagonal.  First,  let 
-1  /2 

V  ■  X  to  diagonalize  the  noise  covariance  matrix.  Then, 


MV=E  M^1/2  X  X+  M^1/2  =  I 


S  =  M‘1/2  S 

V  X  X 


Next,  rotate  the  coordinates  by  a  unitary  transformation  U  so  that  the  S 
vector  is  non-zero  only  in  the  first  component,  and  normalize  its  amplitude 


to  unity 


Y  =  (sf  M"1  S  )'1/2  U  V 
vxx  x' 


*  r;1/2|j  «;,/2  * 


where  rQ  is  the  S/N  ratio  with  optimum  weights,  given  by  (4).  Then, 


"y  ■  u  h;’/2  mx  m;1/2  -  7-  1 

0 


S  =  r_1/2  U  M-1^2  S  = 
y  o  x  x 


Note  that  the  output  S/N  ratio  in  the  new  coordinate  system  is 


sy  <  sy  ■  <Vn 


i.e.,  the  (1,1)  element  of  the  matrix  M‘^. 

The  sample  covariance  matrix  algorithm  can  be  analyzed  conveniently 
in  the  new  coordinate  system.  The  subscript  y  will  be  dropped  in  the 
following  equations.  Again,  the  sample  covariance  matrix  is 


M  ■  A  £  Y.  Y? 


K  ‘k 


The  weights  based  on  this  estimator  are 


W  =  M"'  S 


When  these  weights  are  tested  on  a  different  set  of  samples  than  those 
used  in  estimating  W,  the  S/N  ratio  r1  is 


=  J_w_sC 

r1  *+  ~ 

1  W  M  W 


The  ratio  of  r^  to  the  S/N  with  optimum  weights  is 


Pi  =  r 


1  =  (S"  M'1  S)2 

o  (S~  M'1  S)(SX  M'1  M  M'1  S) 


The  probability  density  of  this  variable  p,  was  derived  in  [1]. 
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In  some  cases  of  interest,  the  weights  may  be  applied  to  the  same 
set  of  samples  used  in  computing  W.  In  this  case,  the  output  S/N  ratio  is 

r  =  S+  M"1  S  =  (M"1 )  -j  (22) 

The  analysis  of  [1]  can  be  extended  to  obtain  an  expression  for  the  pro¬ 
bability  density  of  r. 

n  31 

The  sample  covariance  matrix  has  a  complex  Wishart  distribution1-  ’  , 

i.e.. 


P(A)  =  -Lj-'jjjy-  exp  j-tr(M_1  A)| 


(23) 


where  I  A)  denotes  the  determinant  of  A,  N  is  the  number  of  elements  in  the 
array,  tr  denotes  the  trace  of  the  matrix,  and  A  =  K  M.  The  constant  I(M) 
is  a  function  of  K,  N,  and  the  covariance  matrix  M.  In  (23),  P(A)  is  the 
joint  probability  density  of  the  elements  of  A,  and  is  restricted  to  those 
matrices  A  which  are  positive  definite.  It  assumes  that  the  underlying 
noise  process  is  complex  Gaussian, 

Consider  the  following  representation  of  the  matrix  A. 


(24) 


where  A^  is  a  scalar,  Ag^  is  a  (N-l)xl  column  vector  equal  to  A^»  and 
A^2  is  a  (N-l)x(N-l)  matrix.  As  in  [1],  A  can  be  factored  as  follows 
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Let 


Dn  =  A11"A12  A22  Al2 
°12  =  A12  =  A21 
°22  =  A22 


(27) 


The  Jacobian  of  the  transformation  from  (A-^,  A-^,  A 22)  to  (D-j -j ,  D22^ 

is  one,  so 


p(Dn’D12’D22>  =  °n  i^22^  I(M)  exp{ 


11+D12  °22  D12  +  tr  °22 


)ro} 


=  P(^|  ^  ^(^12  ’ ^*22 ) 


(28) 


where 


P(Dn)  =  C,  Df,"  expj-^D^I 


(29) 


and  is  a  constant. 


Representing  A”1  in  the  same  form  as  (24), 


-1  _  /A11  A12N 
,A21  A22 


(30) 


it  can  easily  be  shown  that  A  =  k- —  .  Since  A  =  K  M  and,  from  (22),  the 

U11 

output  S/N  ratio  is 


r  =  (M_1 =  K  A11  =  i 
“  U11 


(31) 


Let  p  =  r/r  =  —  .  From  (29) 

0  ro  U11 


C?  K'  <  ^ 

P(p)  =  k-N+2  exp  ['  p  \  ’ 


K-N+l 


(32) 


Normalizing  this  distribution,  C2  =  »  anci 


-K-N+l 


■ ;;  C- 1 1  .’>» 


(K-N) I  p 


(33) 


This  is  the  probability  density  function  for  the  normalized  output  S/N 

A 

ratio,  p  =  r/r  ,  when  the  same  samples  are  used  in  M  for  computing  the 
weights  and  for  testing  the  weights. 

From  (33), 

K 


p  = 


K-N 


(34) 


Note  that  p  >  1.  When  the  same  set  of  samples  is  used  for  computing  the 
adaptive  weights  and  for  testing  these  weights,  the  output  S/N  ratio  is 
greater  than  it  would  be  with  "optimum"  weights,  i.e.,  WQ  =  M-1S.  (This 
statement  assumes  that  a  signal  is  not  present  in  the  sample  set,  and 
signal  is  defined  by  the  steering  vector  S.  It  will  be  shown  in  Section  5 
that  this  algorithm  using  the  same  samples  leads  to  a  signal  cancellation 
effect  when  a  signal  is  present.) 

4.  PROBABILITY  DENSITY  OF  MEAN  OUTPUT  NOISE  POWER 

In  the  analysis  of  the  preceding  section,  the  steering  signal  had 
the  form  Sj  =  (1 ,0,0,0, .. .0) .  This  represents  the  usual  sidelobe  can¬ 
celler  problem  where  the  signal  is  present  in  the  main  channel  and  is 
either  absent  or  negligible  in  the  auxiliary  channels.  For  example, 
the  main  channel  may  be  fed  by  a  high  gain  antenna  directed  toward  the 
desired  signal  source  and  the  auxiliaries  fed  by  near-omnidirectional 
auxiliary  elements.  As  formulated  above,  the  magnitudes  of  the  weights 
applied  to  all  channels  vary  with  the  noise  field. 

A  more  conventional  sidelobe  canceller  implementation,  illustrated 
in  Figure  1,  constrains  the  weight  on  the  main  channel  to  a  constant 

value,  e.g.,  unity.  While  the  performance  of  the  two  sidelobe  canceller 
algorithms  is  equivalent,  the  configuration  of  Figure  1  with  constant 
main  channel  gain  is  usually  more  convenient  to  implement.  With  a 
constant  unit  weight  on  the  main  channel,  and  no  signal  components  in 
the  auxiliary  channels,  the  output  signal  power  is  independent  of  the 
adaptive  weights. 


0 


Auxiliary  Elements 


Figure  1. 


Sidelobe  Canceller  with  Unit 
Weight  on  Main  Channel 
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In  the  configuration  of  Figure  1,  the  output  of  the  main  channel 
is  B  and  the  (N-l )  auxiliary  element  outputs  are  represented  by  the 
column  vector  Y,  where  Yy  =  (Y^  • • »Y^_i ) •  The  canceller  output  is 

Z  =  B  -  WfY  (35) 

The  output  noise  power, 

[zj”2  =  BE*"  -  W+YB*  -BY+W  +  W+MyW  (36) 

My  =  vvf 

is  minimized  when 

W  =  M"1  YB*  (37) 

Since  the  output  signal  power  is  independent  of  the  weights,  the 
probability  distribution  of  the  mean  output  noise  power  can  be  derived 
from  the  distribution  of  output  S/N.  First  consider  the  case  of  [1] 
where  the  adaptive  weights  are  based  on  a  sample  covariance  matrix 
and  the  resulting  weights  are  tested  on  a  different  set  of  samples. 

As  derived  in  [1]  the  normalized  output  S/N  ratio,  given  by  p-j  of  (21) 
has  the  probability  density 

p(pl)  =  (N-2)  HK-N+l ) !  (1_VN  2plK  N  1i°  ^  Pi  <  1  (33) 


A 
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Let  rip  denote  the  normalized  mean  output  noise  power,  i.e., 


„  .  lEilfi—  =  i 

1  - 2  P1 

uiV  1 


(39) 


/s 

where  W  are  weights  based  on  a  sample  covariance  matrix  and  WQ  are 
optimum  weights  of  (37).  From  (38)  and  (39) 

(n  -l)N‘2 

P(V  *  (ITTJT^ k'-n+I  J !  ~W\  1  *  ni*  00  (40) 

ni 

The  mean  output  noise  power  exceeds  that  with  optimum  weights  by  the 
factor 

S1  =  KTN  <4’> 

As  before,  K  denotes  the  number  of  samples  in  the  sample  covariance 
matrix  and  N  is  the  total  number  of  channels.  The  number  of  auxiliary 
channels  in  the  sidelobe  canceller  configuration  of  Figure  1  is  (N-l). 

Similar  expressions  for  the  distribution  and  mean  of  the  normalized 
output  noise  power,  when  the  weights  are  computed  and  tested  on  the 
same  sample  set,  can  be  derived  from  the  results  of  the  preceding 
section.  The  density  function  for  the  output  S/N  ratio,  p,  is  given 
by  (33).  Again 
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n  = 


1 


p 


(42) 


From  (33)  and  (42), 


P(n)  = 


kk-n+i 

( K-N ) ! 


0  <  n 


(43) 


and 

-  _  K-N+1  _  n  N-K  .... 

n  =  — g—  =  (1  -  — )  144) 

Note  that  the  mean  value  of  n  is  less  than  unity.  The  weights  based 
on  a  set  of  K  noise  samples  reduce  the  noise  in  that  sample  set  to  a 
lower  mean  power  level  than  weights  based  on  the  true  covariance  matrix 
of  the  noise  process. 

5.  PROBABILITY  DENSITY  OF  OUTPUT  SIGNAL 

As  observed  in  the  preceding  section,  the  sample  covariance  matrix 
algorithm  yields  weights  which  are  more  effective  in  cancelling  the 
noise  in  the  same  sample  set  than  optimum  weights.  These  weights  are 
optimized  for  the  particular  set  of  samples  and  would  not  perform  as 
well  as  the  theoretically  optimum  weights  (based  on  the  true  covariance 
matrix)  against  other  input  noise  samples.  The  weights  are  optimized 


to  minimize  the  output  noise  power  in  the  same  specific  set  of  samples. 
This  suggests  that  these  weights  will  also  tend  to  cancel  any  signal 
present  in  the  sample  set  when  the  weights  are  computed  from  and  applied 
to  the  same  samples.  The  following  analysis  shows  that  this  is  indeed 
the  case. 

Again,  consider  the  sidelobe  canceller  configuration  of  Figure  1, 
where  the  weight  in  the  main  channel  is  constrained  to  unity  and  a 
signal,  if  present,  occurs  only  in  the' main  channel.  Let 

B  =  Yq  +  S  (45) 

where  Yo  is  the  noise  component  and  S  is  the  signal  component  of  B. 

The  weights  based  on  a  set  of  K  input  samples,  (Yk>  in  the  auxiliaries 
and  {B.}  in  the  main  channel,  are 
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The  sidelobe  canceller  output  is  then 

Z  =  B  -  W+Y  (48) 

We  are  interested  in  obtaining  the  probability  density  and  mean 
value  of  the  output  signal  when  a  signal  is  present  in  B  and  the  same 
samples  are  used  for  computing  and  testing  the  weights.  The  mean 
output  signal  voltage  in  Z  is  the  correlation  of  Z  with  S,  i.e., 

SZ*  =  SB*  -  SYfW 

=  SYQ  +  SS*  -  SY-^1  Y(Y*+S*)  (49) 

=  $(Y*  -  Y^M_1YY*)  +  (ST  -  SY^M^YSt*) 

In  (49),  the  first  term  contains  the  product  of  S  and  a  noise  output 
voltage,  so  its  mean  value  is  zero.  The  noise  component  of  this  first 
term  is  the  output  noise  voltage  that  would  be  present  in  the  absence 
of  a  signal.  When  the  signal  is  constant,  the  second  term  reduces  to 


SZ*  =  SS*(1  -  Y+  M"1  Y) 
=  SS*  (1  -  q) 


(50) 


A 
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where  q  is  the  fractional  loss  in  signal.  The  same  expression  of  (50) 
is  obtained  when  the  signal  is  any  function  of  time,  not  necessarily  a 
constant,  provided  its  amplitude  remains  constant.  In  particular,  (50) 
gives  the  mean  output  signal  amplitude  for  a  constant  amplitude  phase 
coded  signal. 

The  probability  density  function  for  q  =  Y  M  Y  is  derived  in  the 

appendix.  The  same  set  of  samples,  { Y^} ,  is  used  in  forming  the  sample 

covariance  matrix  and  in  obtaining  the  mean,  Y.  The  probability 

density  of  q,  assuming  K  samples  of  Y.  are  contained  in  both  T  and  M 

k  y 

and  that  is  a  zero  mean  complex  Gaussian  process,  is 

p(«)  -  rik-wli’rtn-vy in'3(M)k~w  o<q<i  <5i> 

where  N-l  is  the  number  of  auxiliary  elements.  The  mean  value  of  q  is 

q  •  ^  ,  (52) 

N-l 

so  the  mean  signal  voltage  is  reduced  by  the  factor  (1  — ).  This 
is  the  same  factor  obtained  earlier  in  (44)  for  the  reduction  in  mean 
output  noise  power. 

6.  CONCLUSIONS 

In  a  digital  adaptive  array  or  sidelobe  canceller,  the  sample 


covariance  matrix  algorithm  can  be  used  to  obtain  adaptive  weights 
based  on  a  set  of  array  output  samples.  Applying  these  weights  to 
the  same  set  of  samples  to  obtain  the  corresponding  set  of  array 


17 

outputs  is  an  effective  method  of  optimizing  the  system  performance  in 
a  rapidly  changing  noise  field.  The  weights  are  always  optimized  for 
the  specific  noise  field  which  is  present  during  the  sampling  period. 
The  probability  density  functions  for  the  output  noise  power  and  output 
signal  amplitude  have  been  obtained  for  this  same  sample  algorithm.  It 
was  shown  that  the  application  of  same  sample  weights  reduces  the  mean 
output  noise  power  (Eq.  44)  and  mean  signal  amplitude  C Eqs .  50  and  52) 
by  the  same  factor. 


/ 
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APPENDIX 

Distribution  of  X%*TX 

To  obtain  the  probability  distribution  of  this  quadratic  form 
consisting  of  complex  variables,  a  corresponding  quadratic  form  will 
first  be  defined  in  terms  of  real  variables  and  real  matrix  elements. 
Let 


X  =  U  +  i  V 


Z 


(Al) 


where  U  and  V  are  column  vectors  of  real  variables.  The  corresponding 
covariance  matrices  are 


H'- 


(A2) 


Expanding  these  matrices  in  terms  of  U  and  V 

M  =  (U+iV)(U-iV)t  =  Ul7  +  VVT  +  i  (77  -  UVT ) 


=  (A  +  iB) 


A  -B 


(A3) 


1 

"  7 


B  A 


A-2 


where  UU  =  VV  =  A/2 


VU'  =  UV'  =  8/2 


These  relationships  between  the  second  moments  of  U  and  V  follow  from 
the  usual  assumptions  of  zero  mean  normal  distributions  for  the  com¬ 
ponents  of  U  and  V,  with  U  U  =  VjT  and  OL  =  -ujT  . 

mnmn  nm  mn 

We  will  be  interested  in  sample  covariance  matrix  estimates  of  M 
and  H*.  From  the  equalities  of  (A3),  better  estimates  of  A  and  B  for 
the  real  covariance  matrix  can  be  obtained  by  using 


A  =  UUT+  VVT 

B  =  VUT-  UVT  (A4) 


where  K  complex  samples  in  the  estimates  of  M  correspond  to  2K  real 
samples  in  the  estimate  of  H. 

Next  consider  the  inverses  of  the  covariance  matrices  M  and  H. 
Because  of  the  form  of  H,  H~^  can  be  expressed  as 


where 


(A5) 
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CA  -  DB  =  I 
DA  +  CB  =  0 


(A6) 


In  terms  of  the  real  matrices  C  and  D,  the  inverse  of  the  complex  matrix 
M  is 

M"1  =  (C  +  iD)  (A7 ) 

since  M-1M  =  (C+iD)(A+iB)  =  (CA-DB)  +  i (DA+CB)  =  I  (A8) 

from  (A6) 

Next,  the  quadratic  form  ' 


=  UTCU  +  VTCV  -  UTDV  .  (A9) 

The  corresponding  quadratic  form  of  complex  variables  is 
xVx  =  (U-iV)T(C+iO)  (U+iV) 

=  (UTCU  +  VTCV  +  VTDU  -  U^DV)  (A10) 


Since  the  covariance  matrix  M  and  its  inverse  are  Hermitian  ,  the 
imaginary  part  of  (A10)  is  zero. 
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Having  derived  the  relationship  (A10)  between  the  real  quadratic 

form  Z^H~^Z  and  the  corresponding  complex  quadratic  form  X+M“^X,  results 

for  the  real  variables  case  can  be  used  to  obtain  the  probability  density 
~ T  _] — 

and  mean  of  X  M  X  in  (50).  Consider  a  set  of  K  independent  samples 
of  the  variables  X  =  U  +  i V »  where  U  and  V  are  column  vectors  of  N  real 
variables,  normally  distributed  with  zero  mean.  Let 


with  the  sample  mean 

K 

7-r  I  h- 

k=l 


(All) 


and  moment  matrix 


K 

Also,  let  fi,  =  l  ^  (zk-I)(Zk'I)T 
k=l 


(A12) 


=  H  •  Z  ZT  (A13) 

The  elements  of  H|  are  Wishart  distributed  [Eq.  29.9.6,  p.  405  of  Ref.  4], 

A 

and  independent  of  Z.  The  sample  means  are  normally  distributed  with 


zero  mean  and  a  moment  matrix  H'  of  (A2). 
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The  quadratic  form  of  interest,  XTM  X,  can  now  be  related  to  the 
Hotelling  distribution.  Let 


p  ~ T  — 

r  =  (K-l )Z  H1  z 


(A14) 


The  probability  density  function  of  T  [Eq.  29.11.4,  page  409  of  Ref.  4] 


P(T)  = 


2r(|) 


(K-i)N/2r(M)r(|; 


2  -K/2 

+b\]  » T>0 


(A15) 


where  K  is  the  number  of  samples  in  H-j  and  Z,  H-j  is  a  N  x  N  sample 
dispersion  matrix,  and  Z  is  a  column  vector  of  N  sample  means. 


Next,  from  the  matrix  inversion  identity  and  (A13), 

~  -A  'i. 

H_1  =  (H,+Z  ZT)_1=  fi*1  .  J - - - ] — 

1  +  7  H-1  I 


(A16) 


The  quadratic  form  of  interest  is 

-4^-1.  ~Y  A-l  A 
Q  =  X'  M  X  =  l  H  Z 


(A1 7 ) 


In  terms  of  the  sample  dispersion  matrix  H-| 


-t/-.,  h;’z  zt  5 1  V 

q  =  zV  --kr-/z 

1  +2  H.  Z 


Fh-1! 

(i+z^h:1!) 


where  the  probability  density  function  of  T  is  given  in  (A15). 


To  obtain  the  corresponding  probability  density  function  for  Q, 
note  that 


dT  _  (K-l)1/2 

dQ  "  "  2Q1/2(1-Q)3/2 


(A20) 


From  (A15),  (A19),  and  (A20) 


P(Q)  = 


2r(|)(K-i)  2  0-Q)k/2  n  ~- 
(K-1)N/2  r(— )r(|)  ^ 


/ioT 

2Q1/2(1-Q)3/2 


r(|)(Q)N/2'1(i-Q)K/2‘N/2_1 

r(N/2)r(^) 


,  0<Q< 1 


(A21) 


where  the  variables  N  and  K  represent  the  size  of  the  real  matrix  and 
number  of  real  samples  in  H,  respectively. 
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As  noted  in  (A4),  the  effective  number  of  samples  in  the  sample 
covariance  matrix  H  is  2K.  Also,  when  (N-l)  represents  the  size  of  the 
complex  covariance  matrix  and  sample  mean  vector,  the  corresponding 
dimension  for  the  real  variable  case  is  2 (N-l ) .  Replacing  K  by  2K  and 
N  by  2 (N-l )  in  (A21),  the  corresponding  probability  density  function  for 
the  q  of  (50)  is 

P(q)  =  ^-9 - -  0<q<l  (A22) 

r(N-i)r(K-N+i) 


This  density  function  can  be  derived  in  a  different  way  which  does 
not  require  the  analogy  to  the  real  variable  problem.  As  before,  let 


and 


where 


(A23) 


(A24) 


Again,  the  are  column  vectors  of  N  zero  mean  complex  Gaussian 

/\ 

variables.  The  variables  (X(<-X’)  are  also  zero  mean  Gaussian  variables, 
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and  the  sample  matrix  has  the  complex  Wishart  distribution  of  (23) 
with  K  replaced  by  K-l. 

Next,  consider  the  quantity 


R  = 


xV’x 

/S  /N 

xtm1 1 X 


(A25) 


For  any  column  vector  X,  there  is  a  unitary  transformation  on  X  which 
leaves  the  value  of  each  quadratic  form  in  (A25)  unchanged,  but  reduces 

A 

X  to  the  form  (a, 0,0,0. . . ,o)y.  Without  loss  of  generality,  it  can  be 
assumed  that  M=I.  Then 

R  ■  — 4 -  ,  (A26) 

<S1  >11 

and  from  the  analysis  of  Section  3  and  (29), 

rk-1-n  -r 

p(R)  =  fj^rnry r e  *  R>0  (A27) 


The  numerator  of  (A25), 


a  =  K  X+M_1X  =  KX1"  X 


(A28) 


is  the  sum  of  2N  real  Gaussian  variables  and  has  the  probability  density 


,N-1  -a 

Pla>  ■  Tinfr 


a>0 


(A29) 


The  sample  means  are  independent  of  the  sample  matrix  ,  so  R  and  a 
are  independent.  Let 


/N  1  A 

Q  =  x+Mj  X  =  | 


(A30) 


From  (A27)  and  (A30) 


PCQ)  =  Jp(R)P(Q|R)dR 


-J 


RK~N~1  e-R  (RQ)N-1  e~RQ  R  dR 
(K-N-1J!  6  (N-l ) ! R  dR 


(A31 ) 


(K 


(K-l)i 

-N-l)! (N-l)! 


,N-1 


(1+Q)' 


Q>0 


We  are  interested  in  obtaining  the  distribution  of 

^  *1  A 

q  =  X  M  X,  where 


A  A  A  ^ 

M  =  M1  +  X  X1 


(A32) 


Since 


~_1  j.  ~  -1 

.  ,  .  M,  X  X'  M, 

m"1  =  m:1  -  - x-i- 

1+XtM1  *X 


(A33) 


The  quadratic  forms  q  and  Q  satisfy 


1,1  ^ . """ 
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From  (A31 )  and  (A34) 


P(q) 


-(K-.l)  1 


(K-N-l ) ! (N-1 ) ! 


0<q<l 


(A35) 


Replacing  N  with  N-1  in  (A35),  where  (N-1)  is  the  number  of  auxiliary 
elements  as  in  (A22),  the  two  results  for  P(q)  are  the  same. 


